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Abstract 

Heavy-quark interactions in the Quark-Gluon Plasma are analyzed in terms of a selfconsistent 
Brueckner scheme using a thermodynamic T-matrix based on a potential model. The interrela- 
tions between quarkonium correlators, spectral functions and zero-modes, and open heavy-flavor 
transport and susceptibilities are elaborated. Independent constraints from thermal lattice QCD 
can be used to improve predictions for heavy-quark phenomenology in heavy-ion collisions. 
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1. Introduction 

A basic challenge in many-body physics is the understanding of matter properties in 
terms of the forces between the constituents. Medium modifications of the force (or po- 
tential) render its determination an additional challenge. It is therefore important to 
have a good control over the force at least in the vacuum. In Quantum Chromodynam- 
ics (QCD), the fundamental force between static charges, i.e., a heavy quark (Q) and 
antiquark (Q), is well established, both theoretically and phcnomenologically. The heavy- 
quark (HQ) potential has been extracted with high accuracy from lattice-QCD (1QCD) 
computations [1] and can be well represented by the so-called Cornell potential, 

VQQ{r) = -~ + °r, (1) 

characterized by a perturbative color-Coulomb interaction at small distances and a 
"string" term dominant at large r, cf. left panel of Fig. 1. This potential successfully 
describes charmonium and bottomonium spectroscopy in vacuum, which can be under- 
stood as an effective field theory (EFT) of QCD in a I/thq expansion (ttiq: HQ mass). 
The string tension, a ~ 1 GeV/fm, is of nonperturbative origin and most likely associ- 
ated with the gluon-condensate structure of the QCD vacuum. The string term plays an 
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important role already at rather small distances; e.g., for V(r o )=0, i.e., at ro— \ fm, it is 
equal in magnitude (but opposite in sign) to the Coulomb term. Consequently, the char- 
monium spectrum is largely governed by the nonperturbative force (e.g., the ground-state 
binding, Eg^~ 0.6 GeV, collapses to ^0.05 GeV if the string term is switched off). With 
a "calibrated" strong force in vacuum at hand one can study its medium modifications 
and infer from these information about the phase structure and transport properties of 
QCD matter. The analysis of quarkonium production and heavy-flavor spectra in ultra- 
relativistic heavy- ion collisions (URHICs) is aimed at precisely at these objectives, i.e., to 
identify signals of deconfmement and to extract heavy-quark diffusion coefficients from 
the produced medium (see, e.g., Refs. [3-6] for recent reviews). 

The HQ free energy has also been computed at finite temperature, cf. right panel 
of Fig. 1. One observes a gradual penetration of medium effects to smaller distances, 
naturally interpreted as a decrease of the "Debye" screening length, ro ~ l/m/j {mo- 
Dcbye mass). However, even at temperatures as high as 2T C «350MeV, the free energy 
still levels off at a positive value indicative for nonperturbative effects (string term). 
The applicability of the potential approach requires the 4-momentum transfer to be 
dominantly spacelike, i.e., qo ~ q 2 /2uiq <C |<f|. In the vacuum this is satisfied by the 
smallness of the quarkonium binding energy. In the medium the latter is expected to 
decrease further. The thermal momentum scale of a single heavy quark also remains 
parametrically large not too far above T c , p 2 h ~ IvtiqT 3> T 2 (T 2 : momentum transfer 
from the medium). Thus, thermal Q-Q and Q-medium interactions are essentially static 
and elastic but involve nonperturbative interactions. This suggests the possibility of a 
unified description of heavy quarkonia and heavy- flavor transport, with a simultaneous 
treatment of bound and scattering states including resummations. The thermodynamic 
T-matrix approach, which is based on potential interactions, is such a framework, and has 
been successfully applied to electromagnetic plasmas [7]. The in-medium QCD interaction 
is, of course, much more involved than in QED, but the idea of using input from 1QCD 
has revived the potential approach in recent years [8,3]. In this paper we will elaborate on 
the T-matrix approach for open and hidden heavy flavor in the QGP (Sec. 2), and how 
numerical results for spectral and correlation functions, as well as transport coefficients, 
can be tested by thermal 1QCD (Sec. 3). We conclude in Sec. 4. 
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2. One- and Two-Body Correlations in the QGP 



The commonly studied quantity in thermal lattice QCD characterizing the propagation 
of a hadronic current with quantum numbers a is the imaginary-time (r) correlation 
function which is given by a thermal expectation value as 



G a (r,r) = ((j a (r,f)ji(0,0))) 





(2) 



The second equality is a diagrammatic representation for a meson state in terms of its 
free quark- antiquark loop and a 2-body interaction term. The physical information on 
the excitation spectrum in that channel is given by the spectral function in momentum 
space, p a (E,p) = —2 Im G^(E,p). It is related to the euclidean correlator, Eq. (2), via 

oo 

r ( t\ f dE (T? ^ cosh[£(T-l/2T)] 

G a (r,p;T) = J —p a (E,p;T) sinh[E/2T] , (3) 
o 

which illustrates the difficulty in extracting spectral functions from 1QCD "data" of the 
euclidean correlator, since the latter is only obtained for a finite number of r-points on 
a finite interval, 0<t<^. However, using model calculations for the spectral function a 
straightforward comparison to 1QCD data can be performed and constraints evaluated. 
Note that the energy integration encompasses both bound and continuum states. Within 
a HQ potential framework, the T- matrix formalism is thus an ideal choice [9-11]. 

Let us first focus on the quarkonium (QQ) sector. In ladder approximation, and after 
partial- wave expansion, the in- medium T- matrix takes the form [10,11] 

oo 

T a (E;q',q)=V a (q',q) + ^ j dk k 2 V a (q' , k) G° qQ (E; k)T a (E; k,q) , (4) 

o 

where V a denotes the momentum-space potential and Gqq the uncorrelated 2-particle 
propagator (the first diagram in Eq. (2)); its imaginary part can be expressed via the 
in- medium single quark and anti-quark spectral functions, pq and pq, as [12] 

Im G° QQ (E, k) = -J^ (p Q (u, k)p Q (E - co, k)[l - f Q (uj) - f®(E - u)] 

+p Q (oj,k) PQ (E + u,k)lfQ(uj,k) - fQ(E + tJ,k)\) . (5) 

The first term characterizes the standard uncorrelated QQ propagation, while the second 
term (arising from negative-energy contributions) represents Q — > Q (or Q — > Q) scat- 
tering which is nothing but the widely discussed zero-mode contribution to quarkonium 
correlators [13,14]. The scattering effect is encoded in the medium modifications of the 
single-quark spectral function, and thus intimately related to HQ transport. In the vector 
channel (a=V), the pv = 00 component represents the density-density correlation (jo: 
density) which in the static limit gives the HQ-number susceptibility, 

_ d 2 n _dn c _l JdE 2 

Xc(T) " ~W C ~dp~~Tj 2^l-cM-E/T) Pv {E) ' (6) 

o 
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Fig. 2. Selfconsistent Brueckner problem for heavy quarks in the QGP: the heavy-light quark T-matrix 
(upper panel) figures into the calculation of the HQ selfenergy (middle panel) and propagator (lower 
panel), which in turn enters into the 2-particle propagator of the Tmatrix. 

which is also governed by the zero mode. 

The HQ spectral function is determined by the HQ selfenergy in the QGP, Sq, as 

-1 „ 1 



pQ{u),k) 



-Im 



UQ(k) u) — wg(fc) — Sq(u;, k) 



(7) 



cf. the lower panel of Fig. 2. The selfenergy receives contributions from the interactions 
of the heavy quark with the heat-bath particles, £q ~ J f q T q Q, but also from possible 
condensate remnants, especially gluon condensates (recall that the string tension, a, 
seems to survive above T c ), as illustrated in the middle panel of Fig. 2. Following the 
arguments given in the Introduction, the T-matrix approach may also apply to heavy- 
light interactions. In this case one can establish a selfconsistent Brueckner scheme in 
which quarkonium and open heavy-flavor interactions follow from the same potential 
(the construction of such a potential will be discussed below) . Heavy-quark rescattering 
in the QGP is, of course, at the origin of its transport properties. The large HQ mass 
implies that momentum transfers are small compared to its thermal momenta, q <C pth- 
This leads to a Brownian motion picture with a Fokker-Planck equation for the time 
evolution of the HQ distribution function, schematically given by 



d£ 
dt 



7- 



d(pf) , 



dp 



D- 



Qp2 



IP 



f |T 09 | 2 (l-cos( 



(8) 



The thermal relaxation rate, 7, is closely related to the scattering rate Tq = — 2ImEQ; 
while the latter is computed with the forward scattering amplitude, the former includes 
a weight with the scattering angle, 9, signifying the effect of isotropization. 

Let us briefly discuss the construction of the interaction potential, V, in Eq. (4). Orig- 
inally it was hoped to extract it model-independently from 1QCD measurements of the 
HQ free energy at finite T (right panel of Fig. 1), in analogy to the vacuum case (left 
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Fig. 3. Color-average QQ free energy at finite T (left panel) as computed in Nf =2 + 1 1QCD (symbols) 
and pertinent fits [11] using the "Coulomb+string" ansatz, Eq. (9). The T-dependence of the fit param- 
eters is displayed in the right panel for 2 different 1QCD inputs [16,17]. In both cases, the T-dependence 
of a s ~ 0.28 — 0.32 turns out to be weak. (Figures reproduced with permission from [11].) 

panel of Fig. 1). However, at finite T, the entropy contribution in F = U — TS causes 
an ambiguity as to whether the free (F) or internal energy (17) should be used. This 
problem is related to the interplay of time scales for the HQ interaction and for the 
thermal relaxation of the medium, with U and F as limiting cases. On the other hand, 
we note that both F and U are computed from thermal expectation values, as differences 
between a system with and without an embedded QQ pair. A more rigorous approach 
should therefore start from (an ansatz for) a "bare" potential (figuring into the T-matrix 
equation) and then calculate the free (and internal) energies using suitable many-body 
techniques. The resulting "zero-point" functions (F, U) can then be compared to 1QCD 
data and the input 4-point function (potential) tuned for optimal agreement. The thus 
obtained potential can be employed to calculate heavy-quark and -quarkonium properties 
in the QGP as described above. 

The free and internal energies are believed to represent limiting cases bracketing the 
uncertainty in potential models. Instead of a functional parameterization of the lattice 
results a field-theoretic approach accounting for the two main components of the po- 
tential [15] has been proven very useful. With an ansatz for the effective Coulomb + 
confining propagators (including couplings), 



D 00 (k;T) 



1)17 



k 2 + : 



(fc 2 



(9) 



four parameters characterizing the respective interaction strength (a s , tuq) and in- 
medium screening (m_D, rhp) can be adjusted to reproduce the finite- T color-average 
free energy from 1QCD fairly well (see Fig. 3). To improve the applicability in the scat- 
tering regime a Breit correction as known from electrodynamics [18] has been introduced, 
Vboui ~ > Vcoui(l — v\ ■ V2), accounting for the magnetic current-current interaction. This 
modification renders the Coulomb potential "minimally" Poincare invariant [19] and 
the high-energy limit of the T-matrix consistent with the tree-level perturbativc QCD 
(pQCD) result [11]. Note that this extension only applies to a Lorentz-vector interaction 
(Coulomb term), not to the scalar confining term. It thus requires the explicit decompo- 
sition as given by the field-theoretic ansatz, Eq. (9). Further systematic improvements, 
e.g., scrutinizing the uncertainty of retardation effects, need yet to be performed. Com- 
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Brueckner Theory of Heavy Quarks in QGP 
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Fig. 4. Flow chart for a many-body approach to describe properties to heavy quarks and quarkonia in 
the QGP with applications to 1QCD and heavy-ion data. 

parisons to EFT approaches at finite T [20] could also prove useful, even though they are 
usually based on perturbative scale hierarchies which are problematic for the confining 
term. Interesting results are also emerging from lattice simulations of classical Yang-Mills 
theory [21], supporting the importance of nonperturbative effects. 



3. Brueckner Theory of Heavy Flavor in QGP 

In Fig. 4 we display a flow chart envisioning a possible implementation of Brueckner 
theory for heavy quarks in the QGP. Starting from an ansatz for the bare Q-Q potential 
(with suitable corrections and extensions to the heavy-light sector), two-body scatter- 
ing amplitudes are readily calculated. From the heavy-light T-matrix one obtains a HQ 
selfenergy, which has to be iterated for selfconsistency for use in the final T-matrices. 
On the one hand, the latter are employed to compute spectral and correlation functions, 
susceptibilities and free/internal energies (not shown), all of which can be constrained by 
lattice "data" (and used to narrow down the input potential). On the other hand, one 
calculates HQ transport coefficients and quarkonium properties (masses, lifetimes, disso- 
ciation temperatures) which are readily implemented into phenomenological descriptions 
of URHICs (e.g., Langevin simulations or rate equations in a hydrodynamically evolving 
bulk medium) and checked against experiment [22-27]. 

As an example we show in Figs. 5 and 6 results from a recent selfconsistent Brueckner 
calculation including off-shell width effects in the HQ propagators [12]. When using the 
1QCD internal energy as input potential, one finds HQ widths of considerable magnitude, 
reaching up to Fq~0.1-0.2 GeV for on-shell charm quarks. When implemented into the 
charmonium T-matrix, the ground-state (r) c , J/ip) melting temperature (estimated from 
the disappearance of the bound-state peak) is found to be around 1.5 T c , significantly 
smaller than in calculations neglecting width effects (~2T C ). The resulting euclidean 
correlators, including zero modes according to Eq. (5), show fair agreement with 1QCD 
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Fig. 5. Euclidean charmonium correlators (G a ), normalized to the ones with vacuum spectal function 
but with finite-T kernel, R a (r) = G a (r)/G™ c * (r). The T-matrix results at T=1.2, 1.5 and 2T C (dashed, 
dotted and double-dashed curves) [12] are compared to 2-flavor 1QCD computations [28] in the pseu- 
doscalar (r] c ) channel (left; no zero-mode) and scalar (Xco) channel (right; including zero modes). 
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Fig. 6. Left: charm-quark susceptibility (normalized to twice the free massless limit) from selfconsis- 
tent Brueckner calculations with U or F potentials (colored symbols), 1QCD (squares, circles) [29] and 
non-interacting c-quarks with masses m c =l-2 GeV in steps of 0.1 GeV (dashed lines, top down). Right: 
spatial diffusion coefficient for c quarks from the T-matrix with U (pink [22] and red bands [11]) and F 
potentials (green band) [11], LO pQCD (solid line) and AdS/CFT matched to QCD [30] (lower band). 



data (see Fig. 5), even though quantitative improvements are certainly warranted. Still 
within the same framework, the HQ susceptibility, Eq. (6), has been computed and 
compared to 1QCD data (left panel of Fig. 6). For this quantity the width effects lead 
to an appreciable increase over the zero- width limit, providing better overlap with 1QCD 
data. The results of a calculation with the free energy as potential generate a larger \ c , 
since for Fqq the in-medium HQ mass correction, Atuq — Fqq^-^-oo), is about 0.3 GeV 
smaller than for Uqq. Finally, the resulting spatial HQ diffusion coefficients, computed 
from the selfconsistent heavy-light T-matrix within a Fokker-Planck equation (8) with 
D s = T/(mQj), are displayed in Fig. 6. When using the [/-potential, D s turns out to be 
a factor of 3-5 smaller (indicating stronger coupling) compared to pQCD. Toward T c , the 
values are not very far from the strong-coupling limit represented by AdS/CFT (close 
to T c the T-matrix is not yet reliable, e.g., due to the lack of coupled channels such as 
DD). With F as potential D s is roughly in between the [/-scenario and pQCD. 
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4. Conclusions 



We have argued that nonperturbative elastic interactions are key to understanding 
low-momentum interactions of heavy quarks in the QGP. A potential-based T-matrix 
approach is, in principle, capable of connecting heavy quarkonium and heavy-flavor 
transport physics. This opens a rich arsenal of constraints from thermal lattice QCD 
on quantities which directly relate to phenomenological applications in URHICs. Many 
open question remain, e.g., a proper potential definition or the role of correlations near 
T c and nonperturbative Q-gluon interactions. Systematically addressing these issies will 
be essential for a full exploitation of upcoming high-precision heavy-flavor measurements 
at RHIC, LHC and FAIR. 
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